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In order to investigate the effects of nonmagnetic impurities in strongly correlated systems, Quan- 
tum Monte Carlo (QMC) simulations have been carried out for the doped two-dimensional Hubbard 
model with one nonmagnetic impurity. Using a bare impurity potential which is onsite and attrac- 
tive, magnetic and single-particle properties have been calculated. The QMC results show that giant 
oscillations develop in the Knight shift response around the impurity site due to the short-range 
antiferromagnetic correlations. These results are useful for interpreting the NMR data on Li and 
Zn substituted layered cuprates. 



I. INTRODUCTION 

The substitution of impurities into strongly-correlated 
systems is a useful probe of the intrinsic electronic corre- 
lations. For instance, in the layered cuprates it has been 
found that the substitution of nonmagnetic impurities in 
place of Cu in the Cu02 planes strongly influences the 
superconducting [1] and the magnetic correlations [2-7]. 
These experiments have provided valuable insight into 
the interplay of the magnetic, superconducting and the 
density correlations in these materials. In particular, the 
NMR experiments find that when a uniform magnetic 
field is applied, the electronic spins in the Cu02 planes 
in YBa2Cu306+a; polarize to form an oscillatory pattern 
which decays away from the impurities [2-7]. These are 
observed as giant oscillations in the Knight shift response 
of the system. Furthermore, the Knight shifts for nuclei 
close to the impurity have a perfect Curie- Weiss temper- 
ature dependence. The origin of these behaviors is an 
important question. 

Theoretically, the problem of a nonmagnetic impu- 
rity in strongly correlated systems such as the two- 
dimensional (2D) t-J and Hubbard models have been 
studied using various many-body techniques. The single- 
particle properties of an impurity in 2D t-J clusters were 
studied with the exact diagonalization technique [8,9]. 
The magnetic and single-particle correlations around a 
nonmagnetic impurity in the doped 2D Hubbard model 
were studied with the the slave boson and the Quan- 
tum Monte Carlo (QMC) techniques [10]. These cal- 
culations have found Friedel oscillations in the electron 
density. There are also RPA type of calculations of the 
Knight shifts for a nonmagnetic impurity in the 2D Hub- 
bad model [11-13], which suggest that the locally en- 
hanced antiferromagnetic (AF) correlations around the 
impurity are responsible for the giant oscillations and 
the Curie- Weiss behavior observed in the Knight shift 
measurements. 

In this paper, the magnetic and single-particle correla- 
tions around one nonmagnetic impurity in the Hubbard 



lattice will be studied with the QMC technique. For 
this purpose, the determinantal QMC algorithm intro- 
duced by Blankenbecler, Scalapino and Sugar [14] and 
described by Ref. [15] will be used. The emphasis in this 
paper will be on the Knight shift response of the system 
in the metallic state near half-filling. The QMC data 
show that as short-range AF fluctuations grow, giant os- 
cillations develop in the Knight shift response. Here, the 
real-space structure of these oscillations as well as the 
pattern in the electron density will be discussed. In ad- 
dition, the effective impurity potential V e // will be ex- 
tracted from the QMC data on the single-particle Green's 
function. It will be seen that V e ff is weakly attractive 
at the sites neighboring the impurity. Because of the 
"fermion sign problem" , the QMC calculations cannot be 
performed at sufficiently low temperatures where direct 
comparisons with the NMR data can be made. Instead, 
the QMC data will be compared with the RPA-type of 
calculations used for fitting the Knight shift data [12]. It 
will be seen that the real-space patterns of the Knight 
shift, the electron density and the effective impurity in- 
teraction obtained with the QMC technique are in good 
agreement with the results found in the RPA analysis of 
the Knight shift data. These emphasize the role of the AF 
fluctuations in producing the unusual Knight shift data 
on Li and Zn substituted YBa2Cu30g+a;- An important 
feature of these calculations is that here the system is al- 
ways in the paramagnetic state, and no AF droplets have 
formed around the impurity. 

In the following, first the model used for a nonmag- 
netic impurity in the Hubbard model will be discussed. 
Then, the numerical data on the Knight shift and the lo- 
cal electron density will be presented. Next, these results 
will be compared with the previous RPA analysis of the 
Knight shift measurements. Finally, the QMC data on 
the effective impurity potential will be presented. 



1 



II. MODEL 

The 2D Hubbard model is denned by 
H = -t ( c L c j<? + c )a c io) + U J2 ni T n a ~ A 1 n * 

(1) 

where t is the hopping matrix element, U is the onsite 
Coulomb repulsion and \x is the chemical potential. The 
electron creation (annihilation) operator at site i with 
spin a is represented by c\ a (c^), ni„ — c\ a Ci a is the 
electron occupation number with spin a and rij = + 
rin. In the following, numerical data will be presented 
for (n) = 0.875, and U = At and 8t. 

A nonmagnetic impurity located at the origin will be 
modelled by the onsite one-electron potential 

V imp = V a n . (2) 

The value of Vb will be taken to be — 20i, so that it is 
greater than both U and the bandwidth. For this value of 
Vb, the electron occupation number of the impurity site 
is nearly two. In addition, the local magnetic moment 
defined by 

m(r) = V((m*(r))2), (3) 

where m z (ri) = — rijj,, is close to zero at the impurity 
site r = (0, 0). So, for this value of Vb, the impurity site 
is nearly doubly-occupied and magnetically inert. Hence, 
one would expect that this model captures some of the 
physical properties of nonmagnetic impurities in strongly 
correlated systems, even though this form of the bare 
impurity potential is a very simple one. 

The Knight shift of the various nuclear sites is deter- 
mined by the magnetic susceptibility \ through the rela- 
tion 

fc ( r *) = ^2x(r t ,r^iuj m = 0), (4) 

3 

where j sums over the whole lattice [12]. For instance, 
if the hypcrfinc coupling for the nuclear spin Ij at site rj 
is Ali ■ Si, then the corresponding Knight shift is given 
by (7 e /27 n )Afc(r,), where j e and 7„ are the electronic 
and the nuclear gyromagnctic ratios. Here, the staggered 
magnetic susceptibility \ is defined by 

rP 

X (ru rj, iuj m ) = J dr (mT (r)m+(0)) (5) 

where ui m — 2rrnrT is the Bose Matsubara frequency, 
m+ = c| T Cjj, n\ = c^c iT and mj"(r) = e HT m^e~ HT . 
The expectation value (...) is evaluated with respect to 
the hamiltonian of the impure system, which is given by 
H + Vi„ lp . Note that for the translationally invariant 



pure system, fc(r) reduces to the uniform magnetic sus- 
ceptibility Xp«re(q — > 0,io; m = 0). In addition to these, 
results on the local magnetic moment m(r) will be shown. 
Note that the site-dependent electron density (n(r)) and 
m(r) were previously calculated [10]. Here, these quan- 
tities are reproduced for comparisons with fc(r) and the 
effective impurity potential. 

The transferred hyperfine coupling between the nuclear 
spin of a 7 Li impurity, 7 I, and the electronic spins in the 
CuC"2 planes can be described by 

4 

(fl-^Si, (6) 

i=l 

where C is the hyperfine coupling constant and i sums 
over the four nearest-neighbor Cu sites of the 7 Li impu- 
rity. In this case, the Knight shift for 7 Li is 

7 i^i(^)C4fc(r=(l,0)), (7) 

so that the temperature dependence of 7 K is determined 
by k(r) evaluated at the nearest-neighbor site of the im- 
purity, r = (1,0). The measurements of the 7 Li Knight 
shift and of the 89 Y Knight shift at the first and the sec- 
ond nearest-neighbor Y sites of the impurity determine 
fc(r) in the vicinity of the impurity [2,3,5]. The analysis 
of the 7 Li and 89 Y Knight shift data on Li and Zn substi- 
tuted YBa2Cu306+ir concluded that fc(r) develops giant 
oscillations as T decreases [5,12]. The T-dependent line 
broadening of the 63 Cu NMR spectra in Zn substituted 
YBa2Cu306.7 has been also attributed to the develop- 
ment of a staggered polarization of the electronic spins, 
when a uniform magnetic field is applied [7]. The QMC 
data presented here show that it is the AF correlations 
which are responsible for the pattern in fc(r). 

In addition to the magnetic properties, the single- 
particle properties around the impurity will be studied. 
The single-particle Green's function is defined by 

rP 

G(Ti, Tj ,iu> n ) = - dTe^(c ta (T)c] a (0)}, (8) 
Jo 

where u> n = (2n + l)nxT is a Fermi Matsubara frequency 
Here, first the site-dependent electron density (n(r)) will 
be discussed. In order to understand the r dependence of 
(n(r)), the effective impurity potential will be extracted 
from the QMC data on G of the impure system and Gjj of 
the pure Hubbard system. The Dyson's equation relating 
G and Gjj is 

G(r, r', iu n ) = G v {r, r', iu n ) (9) 
+ G u (v,v 1 ,iu: n )T(Y 1 ,Y 2 ,iuj n )G{Y 2 ,T' ' ,iw n ) 

ri,r 2 

where T(r, r', iu> n ) is the effective impurity scattering ma- 
trix. Because of the Coulomb correlations, T(r, r',iu) n ) 
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becomes extended in real space [9]. Here, the diagonal 
component of T is defined as the effective impurity po- 
tential, 

V eff (r,iu n ) = T(r,r,i(j n ). (10) 

In Section III.C, the spatial structure of T(r, r', iu n ) will 
be discussed and results on V e ff (r, mT) will be shown. 

III. NUMERICAL DATA 
A. Knight shift 

The following data were obtained for electron filling 
(n) = 0.875 and an 8 x 8 lattice with periodic boundary 
conditions. First, results will be shown for the U = At 
case at temperatures between 0.25i and l.Ot. In Fig. 1(a), 
fc(r) is plotted as a function of r = |r| in units of the lat- 
tice spacing a. The impurity is located at the origin at 
r = (0,0). For instance, in this notation, r = 1 cor- 
responds to (±1,0) and (0, ±1) sites whereas r = \/2 
corresponds to (±1,±1) and (±1,=f1). In addition, here 
fc(r) is shown in units of t _1 and t is taken to be unity. 

In Fig. 1(a), it is seen that k(r) vanishes at the im- 
purity site. This is because the strong impurity poten- 
tial Vo = — 20t overcomes the electron-electron repul- 
sion U and binds an up-spin and a down-spin electron at 
r = (0, 0), so that the impurity site becomes magnetically 
inert. In the vicinity of the impurity, fc(r) exhibits oscil- 
lations, which grow as the temperature decreases. The 
horizontal long-dashed line in Fig. 1(a) represents the 
value of fc(r) for the pure system at T = 0.25t. Hence, 
away from the impurity, fc(r) goes to its value for the 
pure system, which is the uniform magnetic susceptibil- 
ity. Since these calculations are carried out at high tem- 
peratures, the uniform susceptibility does not yet exhibit 
the temperature-independent Pauli behavior. These data 
on k(r) show that the nonmagnetic impurity strongly in- 
fluences the zero-frequency magnetic correlations around 
it. 

In Figure 1(b), results on m(r) are plotted in the same 
way as in Fig. 1(a). Here, it is seen that, between 
T = 0.25t and l.Ot, m(r) is very weakly dependent on 
the temperature. At r = (0,0), m(r) has a small value. 
Comparing with Fig. 1(a), it is observed that the oscilla- 
tions in fc(r) are stronger than those in m(r). Note that 
both in fc(r) and m(r), the maximum occurs at r = (1, 0). 
The NMR experiments also find that fc(r) has its maxi- 
mum value at r = (1,0). 

In order to understand the origin of the oscillatory 
structure in k(r), it is useful to consider the Fourier trans- 
form of x(r, r',iu) m ) defined by 

X (q, q', iw m ) = ]T e-^ r - q '- r \(r, r', iw m ). (11) 



This quantity is related to the Fourier transform fc(q) = 
J2 r e ~* q r fc( r ) through 

fc(q) = x(q,q' = 0,tu m = 0). (12) 

In Fig. 2(a), — fc(q) is plotted at different temperatures 
as a function of q along various cuts in the Brillouin zone. 
For comparison, in Fig. 2(b) the diagonal susceptibility 

x(q) = x(q,q,^m = o) (13) 

versus q is plotted. In these figures, both fc(q) and x(q) 
are plotted in units of t _1 . Here, it is seen that — fc(q) 
develops a sharp peak at q = (w, n), which closely follows 
the development of the peak in x(q) and, hence, of the 
AF correlations in the system. This relation between 
— fc(q) and the AF correlations becomes more clear if the 
following RPA expression for x(q, q',iw m ) is considered 
[11,12], 

X(q,q',iw ro )= xo(q,q',«w m ) (14) 

+U^2x(<l,<i\iu m )xo(<i\ci',iL; m ), 
q" 

where U is the effective irreducible vertex in the particle- 
hole channel and xo(q",q';^m) is the susceptibility of 
the impure U — system. For q' = and iuo m — 0, one 
obtains 

fc(q) - fc (q) +17^ x(q,q",0)fc (q"), (15) 
q" 

where fc (q) is for the impure U = system. This expres- 
sion shows that the peak in fc(q) at q = (n, tt) is coupled 
to the AF fluctuations. The QMC data seen in Figures 
2(a) and (b) show that this is true for the real system 
also and that it is the AF correlations which cause the 
development of the giant oscillations in fc(r). 

At this point, it is necessary to discuss the physi- 
cal meaning of the q dependence of fc(q). The QMC 
data and the Knight shift experiments show that fc(q) 
peaks at q <~ (tt,tt). Since fc(q) is equivalent to the off- 
diagonal susceptibility x(q, q' = 0,iuj m = 0), this result 
means that the scattering of the AF spin fluctuations 
with ~ (7r, 7r) momentum transfer is one of the impor- 
tant effects of the nonmagnetic impurities. This then 
gives information about the response of the system in 
the magnetic channel to a perturbation in the density 
channel. 

It is useful to compare the structure in fc(r) with the 
site-dependent electron density (n(r)). Figure 3 shows 
(n(r)) versus r for the same parameters as in Figures 
1 and 2. At the impurity site r = (0,0), the electron 
density is nearly equal to two and it is not included in 
this graph. There are oscillations in (n(r)) which decay 
as one moves away from the impurity. For sites close 
to the impurity and especially for r = (1,0) and (1,1), 
(n(r)) increases as T decreases. 
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Figures 4(a) through (d) show QMC data on fc(r), 
— fc(q), x(q) an( i ( n ( r )) f° r the stronger coupling case 
of U = 8t. These data are plotted in the same way as 
for U — At above. The structures in k(r), fe(q), x(q) 
are similar to what have been seen for U = At, only the 
amplitudes differ. When U increases to 8t, the impu- 
rity site remains doubly occupied, however the structure 
in (n(r)) changes. Compared to the U = At case, here 
(n(r = (1,0))) has a bigger enhancement over the back- 
ground. As it will be seen in Section III.C, this has to 
do with the fact that V e ff(r,iirT) at r = (1,0) is more 
attractive for U = 8t. 

These calculations for U — At and 8t were repeated 
for an electron filling of 0.94. In this case, the system 
has longer-range AF fluctuations. Consequently, fc(r) is 
more enhanced and has stronger oscillations. However, 
the structures in fc(r), fc(q), x(q) an d ( n ( r )} are similar to 
what have been shown for the (n) = 0.875 case. Because 
of space limitations these results will not be presented 
here. 

B. Comparison with the RPA analysis of the Knight 
shift data 

The QMC data on the Knight shifts are at high tem- 
peratures to be directly compared with the experimen- 
tal data. However, the real-space structure of k(r) and 
(n(r)) calculated with QMC are in good agreement with 
the results of the RPA analysis of the experimental data 
[12]. In the RPA approach, an effective impurity poten- 
tial which is extended in real space and independent of 
frequency, 

4 

V e ff = V n + V 1 ^2 / n i , (16) 
»=i 

was used to model the nonmagnetic impurity. The onsite 
component Vo was set to a large attractive value, and the 
extended component V\ , which is acting at the four sites 
neighboring the impurity, was used as a fitting parame- 
ter. Both Vo and V\ were taken to be real valued. Then, 
the magnetic susceptibility and fc(r) were calculated by 
using an RPA-like approach for treating the Coulomb 
correlations. The magnetic susceptibility xo of the im- 
pure U = system entering the RPA expression was cal- 
culated by including the self-energy and the vertex cor- 
rections due to the impurity scattering. It was found that 
for an effective quasiparticle bandwidth of leV, the 7 Li 
Knight shift data on the optimally doped YBa2Cu306+z 
can be fitted over a wide temperature range by using 
Vi = -0.15t. 

In this approach, the magnitudes of the quantities such 
as k(r) or the value of V\ used for fitting the data depends 
on the effective bandwidth, and it is difficult to know 
what the exact value of the effective bandwidth should 



be in such a model. In addition, there are uncertainties in 
the values of the hyperfine couplings. For these reasons, 
what should really be compared with is the pattern of 
the oscillations in fc(r) and (n(r)}, and not so much the 
magnitudes. 

In Fig. 5(a), fc(r) obtained by fitting the 7 Li Knight 
shift, 7 K, in optimally doped YBa2Cu306+ir with 7 Li im- 
purities is shown at 100K and 400K. In this calculation, 
the 7 Li hyperfine coupling was taken to be 1.8 x 10 _20 erg, 
corresponding to 0.85kOe//xs. It is useful to compare 
k(r) seen in Fig. 5(a) with the QMC data presented in 
Figures 1(a) and 4(a). Here, it is seen that as T decreases, 
the QMC data develop real-space structure which is sim- 
ilar to the results of the RPA analysis. 

It is difficult to compare the experimental data with 
the QMC calculations directly, but it would still be useful 
to discuss the temperature range where the QMC calcula- 
tions have been performed. The single-band parameters 
which are considered to be appropriate for the cupratcs 
are t ~ 0.45eV and U of order I2t [16]. It is difficult 
to reach low temperatures with large U/t. Hence, here 
calculations have been carried out for U = 8t and At. For 
U = 8t, QMC data have been shown down to T = 0.33t. 
For this value of U, the magnetic exchange J ~ At 2 /U 
is 0.5i, and, hence, T = 0.33i corresponds to two-thirds 
of J, which is too high to make comparisons with the 
experiments. For the intermediate coupling U = At case, 
the lowest T where QMC data are shown is 0.25i. In 
this case, the magnetic exchange J <~ t, and T = 0.25t 
corresponds to J/4. In fact, the QMC data on fc(r) for 
U = At and T = 0.25t, which are shown by the filled cir- 
cles in Fig. 1(a), compare well with the RPA fits to the 
experimental data at 400K, shown by the empty circles 
in Fig. 5(a). This is encouraging, however, the expres- 
sion J ~ At 2 /U for the magnetic exchange is valid in the 
strong coupling limit. So, it is difficult to compare the 
QMC data directly with the experiments. On the other 
hand, studying how fc(r) evolves with U/t and T/t shows 
what to expect for the Knight shifts at low temperatures 
in the strong coupling limit. 

Within the same RPA framework and using an effective 
bandwidth of leV, the Knight shift data on 7 Li and 89 Y 
were also fitted for the underdoped YBa2Cu306+ K with 
nonmagnetic impurities. In this case, it was found that, 
as T decreases from 400K to 100K, k(r = (1, 0)) increases 
from 0.8i _1 to 2.1t~ 1 , whereas k(r = (1,1)) goes from 
0.04i _1 to — 0.8t _1 . This clearly demonstrates the severe 
effect of the nonmagnetic impurity on the zero-frequency 
magnetic correlations in the underdoped cuprates. At 
the temperatures where they are performed, the QMC 
calculations also find that the oscillations in fc(r) are fill- 
ing dependent. For example, for U = At and T = 0.25i, 
fc(r = (1,0)) increases by about 10%, when (n) varies 
from 0.875 to 0.94. On the other hand, if (n) goes from 
0.875 to 0.80, then fc(r = (1,0)) decreases by about 30%. 
In this paper, QMC results are shown for (n) = 0.875, 
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because, for this value of the filling and at the tempera- 
tures where these calculations are carried out, the system 
has short-range AF correlations which are not weak, and 
their effects are visible. 

Next, comparisons are made for the local electron den- 
sity. Fig. 5(b) shows (n(r)) versus r obtained from the 
same RPA analysis at 100K for different values of V\. 
The oscillations seen in this figure are stronger than in 
the QMC data, which were obtained at higher tempera- 
tures. Here, it is seen that as V\ becomes more attractive, 
the electron density at r = (1,0) increases. Comparison 
with the QMC data in Figures 3 and 4(d) implies that 
the effective impurity potential at r = (1,0) is more at- 
tractive for U = 8t than for At. In the next subsection, 
it will be seen that this is indeed the case. 

C. Effective impurity interaction 

In order to understand the spatial structure of (n(r)}, 
it is necessary to discuss the effective scattering matrix 
T(r, r', ioJ n ), which can be obtained from 

T(»W n ) = G^{iw n ) - G-^iiVn). (17) 

Here, T(iui n ) is a matrix of which (i, j)'th element is 
T(ri, Tj, iu n ), T _1 is the inverse of T, and similarly 
for G and Gu- In the following, QMC results will be 
shown for the lowest Matsubara frequency uj n = irT. 
For r and r' more than one or two lattice spacings away 
from the impurity, T(r, r' ,mT) is small in the parame- 
ter regime the QMC calculations were performed. The 
diagonal terms of the scattering matrix T are repre- 
sented by V e ff(r, iw n ). At the impurity site r = (0,0), 
V e //(r,i7rT) is reduced from its bare value of — 20t due 
to the Coulomb correlations, and it also has an imag- 
inary part. For U = At and T = 0.25t, V e ff(r,iirT) 
at r = (0,0) is -18i + iOM, while for U = 8t and 
T = 0.5i, it is — 15t + i2t. Away from the impurity 
site, V e ff(r, iirT) is largest at (±1, 0) and (0, ±1). In ad- 
dition, the effective scattering matrix has large matrix 
elements between the impurity site and the sites neigh- 
boring it. Hence, the presence of the impurity renor- 
malizes the bare hopping matrix elements between the 
impurity and its neighbors. The off-diagonal terms such 
as T(r = (1,0), r' = (0,1), iirT) are small compared to 
V eff (r = (l,0),iirT). 

Fig. 6(a) shows the real-space structure of the real and 
the imaginary parts of V e f / (r, iirT) for U = At and T = 
0.25i, and Fig. 6(b) shows the temperature dependence 
of V e ff(r = (1,0), iirT). Here, it is seen that V e ff at 
r = (1,0) is weakly attractive. The results for U = 8t 
are shown in Figures 7(a) and (b). In this case, V e ff(r — 
(1,0), iirT) has a bigger value. At T = 0.5t, its real and 
imaginary parts are of order —t/A. This is why (n(r = 
(1, 0)) for U = 8t is more enhanced over the background 
than for At. 



In previous analysis of the NMR data [12], an RPA 
type of approximation was used for treating the Coulomb 
correlations and the effective impurity potential was 
taken to be extended in space with the form given by 
Eq. (16), but it was assumed that V e ff is real- valued and 
independent of frequency. Within this model, the fitting 
of the Knight shift data required that V e f / is weakly at- 
tractive at the sites neighboring the impurity. The QMC 
data presented in Figures 6 and 7 show that the real 
part of Veff (r, iirT) at r = (1,0) is attractive, which 
is in agreement with the RPA fitting of the data. The 
QMC calculations also find that T(r, r', iu n ) depends on 
frequency, and, in addition, it has an imaginary compo- 
nent. Since ImT(r, r', ioj n ) is odd in u n , ImT(r, r', iirT) 
should vanish as the temperature goes to zero. In or- 
der to investigate the effects of V e f / having an imaginary 
component, the RPA calculations were repeated using 
Vi(iu n ) = -0.15t (1 + isign(w„)) instead of Vi = -0.152 
in Eq. (16), and keeping all other parameters the same in 
the calculations. It was found that the imaginary compo- 
nent of Vi affects the magnitudes of k(r) and (n(r)), but 
the patterns in these quantities remain similar. These 
comparisons with the QMC data on T(r, r', iuo n ) support 
the form of V e ff, Eq. (16), used in the RPA analysis of 
the Knight shift experiments. 

IV. CONCLUSIONS 

The substitution of nonmagnetic impurities into the 
layered cuprates provides useful information about the 
many-body physics of these materials. In particular, the 
nonmagnetic impurities constitute a perturbation in the 
density channel, and the Knight shift measurements de- 
termine the real-space resolved response of the system to 
this perturbation. In the NMR experiments on 7 Li, 89 Y 
and 63 Cu nuclei in YP^CusOg+x with nonmagnetic im- 
purities, it has been found that, when a uniform magnetic 
field is applied, the electronic spins in Cu02 planes polar- 
ize to form a magnetization with an oscillatory pattern 
[2-7]. In order to understand the origin of this behavior, 
here QMC data have been presented for one nonmag- 
netic impurity in the 2D Hubbard model. These QMC 
data show that, as the temperature is lowered and the 
AF fluctuations grow, giant oscillations develop in the 
Knight shift response around the impurity. The pattern 
of these oscillations and of the electron density were dis- 
cussed here. In addition, the effective impurity potential 
extracted from the QMC data was shown. Since these re- 
sults are restricted to high temperatures, it is not possible 
to compare directly with the NMR experiments. For this 
reason, the QMC data were compared with the results 
of the previous RPA analysis of the NMR data. It was 
found that, as T decreases, the patterns which develop in 
the Knight shift response and the local electron density 
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are in agreement with those found in the RPA analysis 
by fitting the experimental data. 

In spite of these, a number of open questions remain. 
For instance, it is not known whether, within the 2D Hub- 
bard model with only nearest-neighbor hopping, a fit to 
the perfect Curie- Weiss T dependence of the 7 Li Knight 
shift will be obtained at low T between 100K and 400K. 
Perhaps, a second-near-ncighbor hopping or other terms 
are required for describing the pure material. It is also 
possible that the bare impurity potential is more involved 
than the simple form of Eq. (2) used here. Nevertheless, 
the QMC data presented here clearly show the role of the 
AF fluctuations in producing the giant oscillations in the 
Knight shift response. 
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FIG. 1. (a) Knight shift response function fc(r) versus 
r — |r| away from the impurity at various temperatures. Here, 
r is plotted in units of the lattice spacing. These results are for 
U — At and (n) = 0.875 on an 8 x 8 lattice. The long-dashed 
horizontal line denotes the value of fc(r) for the pure system at 
T = 0.25t. (b) Local magnetic moment m(r) = \J {(m z (r)) 2 ) 
versus r for the same parameters as in (a). 
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FIG. 2. (a) Wavevector dependence of the Knight shift re- 
sponse function fc(q) along various cuts in the Brillouin zone, 
(b) Diagonal magnetic susceptibility x(q) = x(q> ^ = 0) 
versus q. These results are for U = 4t, (n) — 0.875 and an 
8x8 lattice at various temperatures. 
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FIG. 4. (a) fc(r) versus r, (b) — &(q) versus q, (c) x(*j) 
versus q and (d) (n(r)) versus r for {7 = 8t, (n(r)) = 0.875 
and an 8 x 8 lattice at various temperatures. In (a), the 
long-dashed horizontal line denotes the value of k(r) for the 
pure system at T — 0.33i. 
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FIG. 5. (a) Knight shift response function k(r) versus r 
at 100K and 400K obtained for the optimally doped YBCO 
with Li impurities. This result was obtained within an RPA 
analysis of the experimental data by using Vi = — 0.15t. Here, 
Vi is the value of the effective impurity potential at the sites 
neighboring the impurity, (b) (n(r)} versus r obtained from 
the same RPA calculations at 100K for different values of Vi. 
In these figures, the long-dashed horizontal lines represent the 
values of k(r) and (n(r)) for the pure system at T — 100K. 
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FIG. 6. (a) Effective impurity potential V e //(r, inT) versus 
r at T = 0.25i. (b) V e// (r, ityT) at r = (1, 0) versus T. Here, 
the filled circles denote the real part of V e ff, and the open 
circles denote the imaginary part. These results are for U — 4t 
and (n) = 0.875. 
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FIG. 7. (a) Effective impurity potential V e ff (r, otT) versus 
r at T — 0.5t. (b) K e// (r,i7rT) at r = (1,0) versus T. Here, 
the filled circles denote the real part of V e //, and the open 
circles denote the imaginary part. These results are for U = 8t 
and (n) = 0.875. 
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